Counting Statistics of an Adiabatic Pump 
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We consider quantum statistics of charge transmitted through a mesoscopic device in the adiabatic 
pumping process. A general formula for the distribution function of the transmitted charge in terms 
of the time-dependent S-matrix is obtained. It is applied to a few simple examples of the pumping 
cycles. We identify coherent pumping strategies which lead to the binomial charge distribution and 
minimize the fluctuations of the pumped current. Conditions for the ideal noiseless quantized pump 
are discussed. 

PACS numbers: 72.I0.Bg, 73.22.-b, 05.45. +b 



c3 



i 

-a 
c 

o 
o 



> 

o 

o 
o 
o 

-I— > 



i 

C 

O 

o 



Adiabatic charge pumping has attracted considerable 
theoretical and experimental interest. It occurs when the 
Hamiltonian of the system is changed periodically with 
time. At the end of the pumping cycle a finite charge 
may be transmitted through the system. Such a charge 
transfer takes place even in the absence of any dc volt- 
age applied to the system. The idea is originally due to 
Thouless [0 , who showed that in certain one-dimensional 
systems the transmitted charge is quantized in the adia- 
batic limit. This mechanism has been proposed for mak- 
ing electric current standards S. 

Recently research has focussed on adiabatic pumping 
through mesoscopic devices |j] M. Most of the work con- 
centrated on the average pumping current and its meso- 
scopic fluctuations || ^| . The general expression for the 
average transmitted charge was derived by Brouwer 0] 
in terms of the time-dependent S'-matrix based on the 
formalism of Ref. Q. The average charge transmitted 
during a cycle depends on the system's conductance and 
is not in general quantized (?|,|^] . The absence of quanti- 
zation is related to the fact that the true Thouless adi- 
abatic conditions can never be achieved for a system 
with infinite leads and hence a vanishing excitation gap. 
The possibility to create electron-hole pairs with arbi- 
trarily small energy leads to dissipation and violates the 
exact quantization ||. Taken together with the discrete 
nature of charge carriers this fact suggests the presence 
of thermal and shot noise in the adiabatic pumps. A 
better understanding of the noise characteristics of such 
devices is certainly necessary if the prospect of their use 
as current standards B is to be taken seriously. 

In this Letter we address the problem of charge count- 
ing statistics in an adiabatic pump. We have greatly ben- 
efited from the extensive work on the current statistics 
of voltage-biased resistors by Levitov et. al. jl(] Q. 

We consider a scatterer connected with the left and 
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right leads, each having n transverse channels 14 1. Such 
a scatterer is characterized by a 2n x 2n unitary scatter- 
ing matrix S, which can be written in the n x n block 
form 



Here r (r') and t (t 1 ) are left (right) reflection and trans- 
mission matrices correspondingly. We shall assume here- 
after that the S'-matrix is a periodic function of time 
r, with the period to, S(t + tq) = S(r). This time- 
dependence is supposed to be adiabatic, i.e. To is much 
greater than the Wigner delay time of the device. We re- 
strict our attention to open systems, i.e. those with the 
dwell time shorter than the inverse mean level spacing in 
the device. Under these conditions one may neglect both 
the energy relaxation and the Coulomb blockade effects 
inside the device. 

The transmission and reflection matrices in Eq. ([l]) can 
be simultaneously diagonalized by the block- diagonal uni- 
tary matrices U(t) and V(r) (see for example Ref. [py) 



S(t) = U{t)S{t)V\i 



(2) 



Here S(t) is a matrix of the form Eq. ([!]) with real diag- 
onal reflection and transmission blocks. The eigenvalues 
of the latter play the role of the instantaneous reflection 
and transmission amplitudes in the channels. The am- 
biguity in the definition of matrices U{t) and V(r) does 
not affect our final results. 

In each cycle a number of electrons Q may pass through 
the scattering region in a direction which depends on the 
detailed form of S(t). The average charge transmitted 
in N pumping cycles was recently given by Brouwer |4J 
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where charge Q is measured in units of the electron 
charge. We have used the representation where the 
transmitted charge is written as half the sum of that 
through the left and through the right leads |15|, thus 
(73 = diag{l, —1} is a Pauli matrix with the n x n block 
structure. Both quantum and thermal noise lead to fluc- 
tuations in the transmitted charge Q. As a result the 
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charge transmitted in one cycle can be described by a 
certain probability distribution. 

Here we obtain the distribution function of transmitted 
charge for an adiabatic pump characterized by a periodic 
time-dependent S-matrix, S(r). More precisely, we cal- 
culate the probability Pn(Q) to transmit the charge Q 
upon completion of N pumping cycles. It is convenient to 
formulate the results in terms of the generating function 
Fjv(A) of the moments of the transmitted charge defined 
as M 



F N (X) = (e iyA ) = / dQP N (Q)e 
For the S'-matrix of the form Eq. (|^) we obtain 



(4) 



det 



l + S x (T)n(T,r')(Si x (r')-Sl(r')) 



where 

S\(t) = exp{-icr 3 A/4}S'(r) exp{i<7 3 A/4} 
The matrix n(r, r') is defined as 

n(r,T') = V^(r)h(T - t')V(t'). 



(5) 
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Here the diagonal matrix h(r — r') is the time Fourier 
transform of n(e) = diag{riL(e), nn(e)}, where nf,rm(e) 
is the energy distribution function of the left (right) lead. 
The operator in the determinant in Eq. (|^) should be 
understood as an operator in the time space as well as 
a matrix in the space of channels. Finally, the integer 
number N defined as 



Nt 

1 fdr f t dU ,dV 
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(8) 



is the contribution to the generating function arising from 
the chiral anomaly (see below). Equation (|^) is the main 
result of this Letter. Before presenting the derivation, we 
shall illustrate its applications on a few examples. 

Expanding lni^(A) to the first power of iX one finds 
for the average transmitted charge, (Q) = J^q QPn(Q), 



(Q) 
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In the absence of an external voltage n(r, t') — ► i/[2tt(t — 
t 1 + ir})} +i/{2%)V^(T)dV{T)/dr for r' ->■ r. Expanding 
S^(t') to the first power in r — r' and using Eq. (^) one 
obtains Brouwer's result Q, Eq. (||). 

A particularly instructive example is given by a single 
channel case (n — 1) with the S'-matrix depending on 
time as 



S(r) = e i ^e(T)/2g e i<T 3 e(T)/2 



(10) 



where 9(t + To) — 0(t) = 27r, and S is time- independent. 
This form of the S-matrix is equivalent to the phase 
winding of the reflection amplitudes, r — > rexp{i9(r)} 
(r' — * r' cxp{— i9(r)}). The contribution from the chiral 
anomaly in Eq. (|8|) coincides with the number of cycles, 
N = N. We concentrate first on a particularly simple 
time-dependence, 9(t) = 2ttt/tq. In this case the deter- 
minant in Eq. (|^) may be easily calculated in the Fourier 
basis leading to 

Fn(X) = J] [l + n L (e k+N )(l - n R (e k ))(e* x - 1)|<| 2 
k 

+ (1 - n L {e k+N ))n R {e k ){e- lX - l)\t\ 2 ]e lNX , (11) 

where e k = n(2k + 1)/(Ntq) are fermionic frequencies. 
At small temperature, T 1/ro, V, this expression sim- 
plifies substantially, leading to 



F N (X) = {e^\r\ z + |t|T(kr + e-^W)~ 



N TQ V 



(12) 



where V is a voltage applied between the left and right 
leads, which is chosen to be such that NtqV/ (2ir) is an 
integer fl0|] . In the absence of external voltage Eqs. (Q) 
and ([12]) lead to the binomial distribution function of the 
charge transmitted through the adiabatic pump 



Pn(Q) = C 



Q N \A 2Q 



(1 



\2\N-Q 



(13) 



Note that only the integer values of the transmitted 
charge have non-zero probability to be detected. The 
physical meaning of this expression is that each pumping 
cycle is associated with an attempt to transfer one elec- 
tron. The success probability of such an attempt is given 
by the reflection probability, |r| 2 , whereas the probabil- 



ity of failure is 1 



The above statistics should be 



compared with the case of the dc voltage applied across 
the scattering region (no pumping). This case may be 
obtained from Eq. (|l^) in the limit TV — > 0, whereas 
NtoV/(2tt) — > N - integer. We immediately recover the 
familiar results fllGl 



Tr{S(T)n(T,T> 3 ,SV)]} +N. (9) 

T' — >T 



PM = C%\tfQ(l-\t\ 2 ) N -Q. (14) 

The distribution is also binomial, however the probabil- 
ity of success is given by the transmission probability, 
\t\ 2 . The second cumulant of the transmitted charge co- 
incide for the adiabatic pump and applied dc voltage and 
is given by 



«Q 2 )H|£| 2 (l-|t| 2 )iV, 



(15) 



leading to the maximal noise power in both cases at 



1^ 



1/2. 



The binomial distribution, Eq. ( pL3| ) was derived above 
for the simplest time dependence of the form r(r) = 



2 



r exp{2niT / tq} . The same result may be shown to be 
valid for a more general class of pumping strategies, 
which we call the coherent pumping, following the ter- 
minology of Ref. 0,0: 



«(t) 



e 27riT/r _ z 

1 — z*e 27TiT / T ° 



(16) 



where z is a complex number with \z\ < 1. Indeed, at 
zero temperature n(efc) = 9(—2k — 1) has infinitely de- 
generate eigenvalues. One can therefore diagonalize the 
operator in Eq. (|^), in the z-dependent basis which does 
not mix positive and negative frequencies. The value of 
the determinant does not depend on z and one recov- 
ers the binomial distribution ([l3]). Such coherent pump- 
ing strategy guarantees the minimal possible value of the 
noise in the transmitted charge. The general expression 
for the noise of the adiabatic pump may be derived by 
expanding lnF^A), Eq. (|^), to the second order in iX 
and is given by 



Nt„ 



«Q 2 )) 



drdr' Tr {l - [S^a^S) r {S^a^S) r , } 



(4Nt ) 



sin 2 [7t (t - t')/(Nt )} 



(17) 



Substituting the S-matrix of the form Eq. (|1C|) and min- 
imizing the charge fluctuations, ((Q 2 )}, with respect to 



[13], one finds that the coherent pumping, Eqs. (10) 



JS(t) 

(pf), lead to the minimal possible noise. The value of 
this minimal noise is given by Eq. (0). Note also that 
the second moment, Eq. (0) (as well as the first one, 
Eq. (^), and the higher ones) may be expressed through 
the S'-matrix only, rather than through auxiliary matri- 
ces defined in Eq. (^). 

Next we consider a 2 x 2 scattering matrix with real re- 
flection and transmission amplitudes given by r = — r' = 
cos(27tt/to) and t = t 1 = sin(27rr/ro) respectively. In 
this case U(t) — V(r) — 1 leading to N = 0, and h is 
the equilibrium distribution function. One can show that 
F/v(A) in Eq. (@) is an even function of A, and is therefore 
real. As a result one may employ the method of Ref. |12j ] 
to compute the determinant of the operator in Eq. flq): 
one multiplies this operator by its Hermitian conjugate 
and takes the square root. The resulting operator may 
be written as 1 + (1 — h)S\S-\n + nS_ x S\(l — h). In 
the energy representation it has a finite number of the 
off-diagonal matrix elements and its determinant can be 
straightforwardly evaluated. This way one obtains 



F N (X) 



1 + cos A 



JV 



(18) 



Therefore such pump is a realization of the "random walk 
motion" for charge. Indeed, there are three possible val- 
ues of the transmitted charge in each cycle: Q = with 



the probability 1/2, and Q — ±1 with the probability 
1/4 each. 

We note that the logarithmic derivative 
iV^ (r)dV (t) / dr is analogous to the instantaneous ma- 
trix of "voltages" applied to the incoming channels. The 
integral of this quantity can be interpreted as the number 
of transmission attempts [0. The different outcomes of 
such attempts lead to the noise of the pumping current. 
In general the probability distribution of the transmitted 
charge is not binomial. If the "voltage" matrix can not 
be diagonalized simultaneously with the reflection and 
transmission matrices then the distribution function of 
the transmitted charge does not factorize into binomial 
distributions of elementary transmission processes. 

In contrast, the matrix U(t) corresponds to the out- 
going channels and enters the final expression d) only 
through the chiral anomaly term ^) and therefore con- 
tributes to the average current but not to the noise. For 
example, the pumping cycle of the form S(r) = U(t)S 
at zero temperature would produce a noiseless quantized 
pumping current. 

We turn now to the derivation of Eq. (@) . To this end 
we model the leads by a 2n-component vector of chiral in- 
coming fermions (jif>r,(x, r), ipii(x, r)) and 2n-component 
vector of chiral outgoing fermions (£l{x, t), £r(x, t)). 
The action for e.g. left lead is written as 



for [dx$ L (dt+v L d x )il>L+ZL(dt-VLd x )ZL , (19) 



where vl is a diagonal nx n matrix of the left lead chan- 
nel velocities. In this expression the time integral runs 
along the Keldysh contour, C, from r = to r = Ntq 
and then back to r = 0. The right lead is described by 
the similar action with the space integral running from 
x = to x = +oo, and the velocity matrix vr. Finally 
the incoming and outgoing channels at x — are related 
by the time-dependent S'-matrix operator 



£(0,r) =v 1 I 2 S{t)v- x ' 2 ^{Q,t). 



(20) 



The current operator has a form / = (Ij,+7r)/2, where 

II (t) = [$ L (0" , t)v l ^ l (0- , r) -£ L (0- , t)vUl (0~ , r)] . 

(21) 

The operator of the charge transmitted in N cycles is 

Nt 

given by Q — J drJ(r). Finally, the generating func- 
o 

tion may be written as 

r -Sl-Sr + % JdT\(r)I(T) 

F N (X)= e c ; (22) 



where A(r) is equal to A on the forward and —A on the 
backward part of the Keldysh contour. The fermion fields 
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in this integral obey the boundary condition, Eq. (20). 
One has to specify the initial, r = 0, density matrix, 
which implicitly defines the Green functions. We fix the 
occupation numbers in the incoming channels of the left 
and right leads to be «i(e) and n^{e) correspondingly, 
whereas the outgoing channels arc supposed to be ini- 
tially empty in accord with the scattering setup. 

The subsequent calculations amount to the evaluation 
of the Gaussian integral in Eq. (|22|). To this end we first 
make the chiral gauge transformation of the fermionic 
fields: ip(x,r) — » V(t)^(x, t) and £(x,t) — » J7(£)£(x,t). 
As a result, the boundary condition for the new fermions 
contains the S(t) matrix only and the action acquires an 
additional time-dependent (matrix) chemical potential 
term 

r dV 



5S= dr Idxij) 







4> + t 



IP 



du 



(23) 



product of binomial distributions. We point out the ma- 
jor role played by the chiral anomaly contribution to the 
average transmitted charge. Such anomalous term did 
not arise in the context of voltage-biased systems Jl0| , 
but is extremely important for the adiabatic pumping 
setup. 
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Such potential term results in the redefinition of 
the density matrix according to Eq. (Q). Impor- 
tantly, upon the chiral gauge transformation the ex- 
pression for the current acquires an extra term / — > 
/ + l/(4iri)Tr{dW(T)a 3 U(T)/dT-dV^T)a 3 V(T)/dT} 
arising from the chiral anomaly |l7j . 

Since the source field, A(r), is a constant on the both 
branches of the Keldysh contour, one may eliminate the 
/ XI term from the action by the time-independent gauge 
transformation [[[o|, e.g. ipL — > e ld ( x ~° ^/ 2 tpL on the for- 
ward branch of the contour and ipL —* e~ l9 ( x ~° '^ 2 ipL 
on the backward branch. Such transformation leads to 
the change in the phase of the forward scattering am- 
plitude and can be taken into account by a redefinition 
of the S-matrix in the boundary condition, Eq. (20), 
S — * S±\ on the forward (backward) branches, with S\ 
defined in Eq. @. 

The subsequent steps are straightforward. One inte- 
grates out all degrees of freedom except for those which 
reside directly at the scatterer, x — 0. Using the bound- 
ary condition, Eq. (20), with S^^matrix, one eliminates 
the incoming degrees of freedom, ij){x — 0,r). The re- 
maining Gaussian integral over the outgoing fermions, 
£{x = 0,t), can be straightforwardly evaluated, result- 
ing in the determinant written in Eq. (|^). The remain- 
ing term, exp{z7VA}, is the contribution from the chiral 
anomaly, as explained above. 

To conclude, we have derived a general expression for 
the counting statistics of the charge transmitted through 
a system described by a time-dependent S'-matrix. The 
only limitations of our result are the requirements of adi- 
abaticity and the absence of inelastic processes in the 
scattering region. The absolute minimum of the noise 
power may be achieved by the coherent pumping strat- 
egy, in which case the charge distribution is given by the 
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